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Abstract 

SPOrt is an ASI-funded experiment specifically designed to measure the sky polar- 
ization at 22, 32 and 90 GHz, which was selected in 1997 by ESA to be flown on 
the International Space Station. Starting in 2006 and for at least 18 months, it will 
be taking direct and simultaneous measurements of the Stokes parameters Q and U 
at 660 sky pixels, with FWHM = 7°. Due to development efforts over the past few 
years, the design specifications have been significantly improved with respect to the 
first proposal. Here we present an up-to-date description of the instrument, which 
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now warrants a pixel sensitivity of 1.7 /iK for the polarization of the cosmic back- 
ground radiation, assuming two years of observations. We discuss SPOrt scientific 
goals in the light of WMAP results, in particular in connection with the emerging 
double-reionization cosmological scenario. 
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1 Introduction 

The last decade of the past century has seen a number of great achievements 
in observational cosmology. Among them, the assessment of the nature of the 
cosmic microwave background radiation (CMBR), by COBE (Mather ct al. 
1990, 1999), and the detection of its tiny temperature fluctuations (Smoot ct 
al. 1992; Bennett et al. 1996), providing a sound confirmation of the gravita- 
tional instability picture. After COBE, a number of ground and balloon exper- 
iments (see, e.g., De Bernardis et al. 2000, Hanany et al. 2000 for a review), 
together with the ongoing WMAP mission (Bennett ct al. 2003), extended 
the measured angular spectrum of temperature fluctuations up to / ~ 800. By 
fitting model predictions to these measurements, the multi-dimensional pa- 
rameter space of available cosmological models was significantly constrained. 

The main target of further analysis was clearly polarization. Besides conveying 
independent information on the cosmological model, so as to remove residual 
degeneracies between various parameters, large angle polarization data allow 
to explore the thermal history of the Universe at low redshift, in the epoch 
when primeval fluctuations gradually turned into today's systems and objects. 
Upper limits to the polarization signal had been put by ground experiments 
(see, e.g., Keating et al. 2001; Hedman et al. 2002), while DASI (Kovac et al. 
2002) claimed a detection with significance of 4.9(7 on angular scales ~ 0.5°. 

The entry into the new century of CMBR exploration was however set by 
WMAP detecting a clear signal of TE correlation at large angular scales 
(Kogut et al. 2003). This was much more than an aid to remove parame- 
ter degeneracies. Its intensity carried unexpectedly important information on 
the history of our Universe, operating a radical selection among reionization 
scenarios. 

This finding makes the SPOrt experiment quite timely. This experiment, fully 

funded by the Italian Space Agency (ASI), was proposed to the European 
Space Agency (ESA) in 1997 to put a polarimeter onboard the International 
Space Station (ISS). Its frequency range (22-90 GHz) and angular resolution 
(FWHM = 7°), together with an almost full-sky coverage, are optimal to 
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search for CMBR polarization on large angular scales. This strategy has now 
received an unexpectedly strong support by the first-year WMAP release, 
which favours values of the Universe optical depth r ~ 0.17. 

Such high values need a confirmation, which shall come first by further WMAP 
releases. An independent confirmation, based on an instrument devoted to 
polarization and operating independently from anisotropy, will then be allowed 
by the SPOrt experiment. It should be stressed that a safe assessment on a so 
large r is not only quantitatively relevant, but changes our views on the history 
of cosmic reionization and on the very epoch when the first objects populated 
the Universe. Besides allowing to measure r, large angle polarization will also 
put new constraints on the nature of Dark Energy (DE), which determines 
the relation between time and scale factor in the epoch when reionization 
occurred. 

Only full sky experiments, based on space, can provide information on large 
angular scale polarization, ground based or balloon experiments being intrin- 
sically unable to access such scales. 

In the light of such considerations, here we report on the present advancement 
of the SPOrt experiment, on the technical aspects of the polarimeter we shall 
fiy and on the physics we expect to explore through its observations. 

The plan of the paper is as follows: Sect. 2 reviews large-scale CMBR polar- 
ization theory as well as the current scenario, including foregrounds; Sect. 3 

describes the SPOrt experiment, with special emphasys on its technological 
and design aspects; Sect. 4 is focused on the main physics results SPOrt is 
expected to provide, compared to WMAP. The conclusion are summarised in 
Sect. 5. 



2 Scientific targets: theory and current scenario 

2.1 Anisotropy and polarization: definitions 

SPOrt aims at measuring the polarization of CMBR. Quite in general, to define 
the properties of linearly polarized radiation coming from a sky direction n, 

we make use of a basis (ei, 62), in the plane perpendicular to n, to define the 
2x2 tensor lij = (EiEj) /Anc {Ei are components of the electric field vector 
E, perpendicular to h). The total intensity of the radiation from the direction 
n is then the tensor trace I = In + l22- By angle averaging, we obtain the 
average radiation intensity /. 
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This average value ought to be subtracted from the signal, to detect the main 
CMBR features. We therefore replace lij by the tensor = Lij/I — 5ij/2. 
Its components are directly related to the temperature anisotropies and the 
Stokes parameters accounting for linear polarization, as follows: 

T ^ (An + A22)/4 , Q = (An - A22)/4 , U = A12/2. (1) 



Q and U depend on the basis (61,62). If such basis is rotated by an angle 
V^, so that e'l = cosV' ci + sin'^ 62 and e'2 = — sin-i/^ ci + cosifj 62, the Stokes 
parameters change into Q' — cos2ip Q + sin2ip U, U' — — sin 2-0 (5 + cos 2-0 U. 
According to Tegmark & de Ohveira-Costa (2001), it is then convenient to 
define a 3-component quantity 

D = (T,g,[/) (2) 

which transforms into 

D' = R(V^)D , (3) 



when the basis is rotated by ip] here 



/I 






cos 2il} sin 2%l) 
— sin 2il) cos 2%l) 



(4) 



is a 3 (8) 3 matrix. 



2.1.1 Angular Power Spectra 

T, Q and U depend on the sky direction n. Prom them we obtain angu- 
lar power spectra (APS), by expanding on suitable spherical harmonics. The 
spherical harmonics currently used to expand anisotropies (T) are however 
inappropriate for polarization, as the combinations Q ±iU are quantities of 
spin ±2 (Golberg et al. 1967). In this case, one makes use of the spin-weighted 
harmonics ±2^™ (Zaldarriaga & Seljak 1997) 

{Q + iU){n) = Eim a+2,im +2yim{n) (5) 

{Q - iU){n) = Y.lma-2,lm -2Ylm{n). 

In order to relate CMBR polarization to its origin, instead of using the Stokes 
parameters Q and U, it is more convenient to describe linear polarization 
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through two suitable scalar quantities E and B (see: Seljak & Zaldarriaga 
1997 and Kamionkowski ct al. 1997a, who give a shghtly different notation). 
Theory provides direct predictions on such quantities and the treatment in 
terms of E and B takes full account of the spin-2 nature of polarization. Their 
scalar harmonics coefficients are related to a±2,/m by the linear combinations 

O-EM — —{<l2,lm + a-2,/m)/2 and aB,lm — {<l2,lm — a-2,/m)/2i 

These combinations transform differently under parity. Thanks to that, just 
four angular spectra Cyi {Y — T, E, B and TE) characterize fluctuations, 
if they are Gaussian. These APS are the transforms of the autocorrelations 
of T, E and B and the cross-correlation between E and T; invariance under 
parity grants the vanishing of the cross-correlations between B and any other 
quantity, and the related APS. More in detail, such APS are defined as follows: 



{(^*B,lm^B,lm') 



-i'Cti] 

-i'CbV, 



^m,m'CTEl- 



(6) 



In order to fit data with models, one must also be able to invert these rela- 
tions, working out signal correlations from APS. This operation can be fairly 
intricate, because of the complex properties of spherical harmonics. One can 
start the procedure by defining the following functions: 



(/-2)!(/-m)! 
\ (Z + 2)!(/ + m)! 



F2,lm{9) — 2^ 



{I - 2)\{l - m)\ m 
\\ (/ + 2)!(/ + m)! 





— I 
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sin^ e 




^sin^^^'- 
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sin^^ ' 



I l{l 



{I 



cos9Pr{cos9) + + m)Pi"l^{cos9)]. (7) 
Then, simple but lengthy calculations provide relations, similar to the equation 

(8) 



(T(l) T(2)) = ^^i±lcTzP/(cosai2) 



47r 



(ai2 is the angle between the directions 1 and 2), relating fluctuations to 
spectra, for Stokes parameters and their correlations. Let us remind that the 
l.h.s. of Eq. (8) prescribes an averaging operation which, in principle, is an 
ensemble average. To work out such function from data, the average is taken 
on all pairs of sky directions with angular separation ai2. 
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The relations analogous to Eq. (8), for the Stokes parameters, arc simpler if 
the basis vector ei is aligned with the line connecting the points 1 and 2. There 
they read 



9/ + 1 

{Qil)Q{2))^J2^^[CEiFi,i2{ai2) - CbiF2,i2M] 

{U{1)U{2))^J2'^^[CbiFi,i2M - CeiF2,i2M] 
I 47r 

(r(l)g(2)) = - ^ ^-^CteiF,,ioM 

(T(l)f/(2))=0. (9) 

For a generic basis, such relations are more complicated (Ng & Liu 1999). Ow- 
ing to their definition and to Eq. (3), however, it is evident that the correlation 
are components of the 3 ® 3 matrix {DiDj), which transforms as follows: 

(dim:) = Rn,nM{DnDj,)R,,j,{^2) (10) 



with R('0) defined in Eq. (4); ■0i_2 are the rotation angles to put the local 
polar basis in 1,2 along the great circle connecting the pair. Making use of 
these transformations, we can pass from the priviledged basis, where Eqs. (9) 
hold, to the generic basis needed for observational or modellization aims. 

When real data are considered, an extra complication cannot be avoided, as 
detected signals suffer from both the finite resolution of antennae and pixila- 
tion effects. The total result is an effective reduction of the actual power at a 
given multipole, and the measured APS become 

^ Cy,i\sWi\' (11) 

with 

sWi^bpi^,isBi (12) 



where bpi^^i is the spherical harmonic decomposition of the pixel window func- 
tion, and is the beam window function. This writes 

sBi^exp[-{l{l + l)-s^)ay2], (13) 



where a = FWHM/v81n2 for Gaussian beams and s = 2 for polarization, 
whereas s = for anisotropy (Ng & Liu 1999). The rms polarization Prms 
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measured over a finite-resolution map is then given by 

PLs-l:^^C,^,Wf^ (14) 

l>2 

where 

Cpi — Cei + Cbi- (15) 

When Q and U maps are available one can also use P^^^ — Q^ms + ^rms)- 

As an example, let us consider an instrument with a FWHM aperture of 7°. 
In this case the reduction caused by a factor {2Wi/bpix,iy is for I ~ 27. It 
is then clear why, with such an aperture, it makes sense to explore spherical 
harmonics up to / = 30. Let us further recall that, in order to obtain an 
optimal exploitation of the signal, it is convenient to set pixels with centres 
at a distance - FWHM/2. 



2.1.2 Small sky areas and non-Gaussian signals 

The treatment of the previous subsection can be simplified, when a small 
patch of the sky is considered, so that the celestial sphere can be treated as 
flat (Seljak 1997). Then the spin-weighted harmonic expansions (5) reduce to 
Fourier expansions in the {9x,0y) plane: 

Yil) = N;^J2^M-^1- 0k)Y{9k) ; (16) 

k 



here A'pix is the number of pixels, and Y stands for either Q or U. Passing to 
E and B is also simpler, as 



B{1) = -Q{f} siTL{(t>i) + U{1) cos{(Pr) (17) 

(pj- being the angle between / and the 9^ axis. From the expansion coefficients, 
the APS can be recovered and read 

Cyi^n{Y{l)Y*{l)). (18) 



Here Q is the size of the sky patch and the average () is to be performed over 

— * 

those I having magnitude close to I. These spectral components, referring to 
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the portion of the sky we consider, are local and can depend on the setting of 
such portion. 

When deahng with foregrounds, which are non-Gaussian, APS can still be 
useful to evaluate their contamination to the CMBR signal. In such a case, 
an angle averaging should be explicitly performed. This amounts to summing 
upon m and dividing by 2/ + 1. Eqs. (9) should not be used, because of the 
phase coherence of harmonic amplitudes. Furthermore, one should carefully 
distinguish between local and global APS. Foreground subtraction ought to be 
performed using local APS; in principle locality should be enhanced as much 
as data uncertainties allow, in order to perform correct subtractions on small 
angular scales. 

It was not recognized for some time that, by expanding the polarized intensity 
n = v^Q^~+TP in spherical harmonics, a further spectrum, Cn/, is constructed. 

As first noted in Tucci et al. (2002) Cjn conveys an information distinct from 
Cpi. In particular, Cui provides an incomplete information, since it does not 
include the polarization direction; for instance, it cannot be used to build sim- 
ulated polarization maps, and cannot account for beam-depolarization effects. 
However, it is still useful when studying the emission of Galactic foregrounds 
at low frequency (< 2 GHz), where Cei and Cbi can be affected by Faraday 
rotation, whereas Cm turns out to be quasi-independent of it (Tucci et al. 
2002), and can be safely extrapolated to higher frequencies. 



2.2 Origins of anisotropics and polarization of the CMBR 

Photons undergoing Thomson scattering become polarized. In the presence 
of uniform and isotropic electron and photon distributions, however, angular 
averaging smears out polarization. To preserve some of it, not only must the 
photon distribution be anisotropic, but anisotropy must have a quadrupole 
component. Until photons and baryons are tightly coupled so as to form a 
single fluid, any quadrupole in the photon distribution is erased; accordingly, 
all scales entering horizon before recombination develop a quadrupole only 
after decoupling begins, so that photon and baryon inhomogeneities are out 
of phase. 

When photons have acquired a mean-free-path of the order of fluctuation 
sizes, therefore, quadrupole anisotropics arise. Such a long mean-free-path, 
however, also means that the optical depth r of baryonic materials is ap- 
proaching unity (from above) and more and more photons are unscattered 
until present. However, only the restricted fraction of photons that still un- 
dergo a late scattering, can acquire polarization. Therefore, primary polariza- 
tion is small both because quadrupole is small and because few photons are 
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scattered after it arises. Quadrupolc (and, hence, polarization) turns however 
higher on scales entering recombination with larger dipole, i.e. in the kinetic 
stage. The number of photons scattered decreases further for scales entering 
the horizon after recombination, where the polarization degree rapidly fades 
with scale increase. 

All these effects can be followed in detail by numerical linear codes, like CMB- 
FAST ^ , able to compute the amount of polarization produced on each scale 
as a function of the parameters of the underlying cosmological model. The 
information that the tiny polarization signal carries can therefore help to re- 
duce any parameter degeneracies that the stronger anisotropy signal cannot 
resolve. In particular, the detailed history of recombination, occurring at red- 
shifts z ~ 10^, can be inspected through these signals. 

However, it is easy to see that the horizon scale at recombination occupies 

an angular scale ~ 10~^ radians, on the celestial sphere. Accordingly, the 
primary polarization spectrum can be expected to peak at / ~ 10^-10^ and to 
be negligible at / < 30, those /'s we can inspect if FWHM ~ 7°. The above 
inspection is therefore precluded to the SPOrt experiment. 

The polarization signal SPOrt aims at detecting must be related to later 
events, when the scales entering the horizon subtended angles ~ 7° or were 
even larger than so. Should the Universe be filled with neutral gas from re- 
combination to now, no such polarization could be expected. It is known, on 
the contrary, that diffuse baryonic materials, today, are almost completely 
ionized. But this could not be enough: baryonic materials are so diluted, to- 
day, that the scattering time for CMBR photons tg ~ 4.45 ■ VS^'^^'^^^h^'^'^ {h: 
Hubble parameter in units of 100 km/s/Mpc) widely exceeds the Hubble time 
til — 3.09 • 10^^/i~^s. Hence, secondary polarization can arise only if reioniza- 
tion is early enough not to be spoiled by baryon dilution. 

Let us therefore extrapolate these figures to the past: while tg oc (a: the scale 
factor), the dependence of tu on a depends on the model. Numerical compu- 
tations are avoided by assuming an expansion regime dominated by ordinary 
non-relativistic matter, as in a standard CDM model. Then an ionized Uni- 
verse is opaque (r = 1) to CMBR a.t z > Zop ^ h {Vtbh)-"^/^. If fi^/i^ < 0.022, 
as from standard BBNS limits (sec Dolgov 2002 for a recent review), and 
h ~ 0.5, we obtain Zop > 64. The redshift range where reionization may have 
occurred (zri) is safely below such Zop. A standard CDM model disagrees with 
data, but no reasonable cosmology allows a substantial reduction of Zop- A full 
rc-scattering of CMBR at low redshift is therefore almost excluded, but the 
opacity r can be non-ncglcgiblc if Zri is not lower than ~ 8-10. In turn, this 
may cause secondary polarization. 



www.physics.nyu.edu/matiasz/CMBFAST/cmbfast.html 
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2.3 Relevance of secondary polarization 



Let us briefly remind which kind of physics is explored through r measure- 
ments and, hopefully, through a more sophisticated analysis, leading to es- 
timates of Zri and other related parameters. They could shed new light on 
(i) the cosmic reionization history and the physics of systems causing it; (ii) 
the expansion laws during reionization and the cosmic substance ruling the 
expansion. 

Good reviews on point (i) were provided by Loeb & Barkana (2001), Shapiro 
(2001), MaseUi et al. (2003) and Haiman (2003). Most studies, assuming the 
flat A-dominated concordance model {flrn = 0.3, h = 0.7), until recently 
converged in setting Zri in the interval 8-10 and suggested that hydrogen 
ionized then and remained such until now. Reionization was attributed to 
radiation produced by metal-free star formation in early galaxies, whose effects 
depend on several parameters, as the photon escape fraction f^sc and the 
slope of the Salpeter Initial Mass Function (IMF). Among the limits that such 
theories must fulfll is that helium, at variance from hydrogen, reionized much 
more recently. However, although leaving all parameters vary in the widest 
possible range, Zri values above 10 had not been envisaged. Other possible 
sources of reionization had been identifled in AGN (Ricotti & Ostriker 2003 
and reference therein) as well as exotic sources like decaying heavy sterile 
neutrinos (see Hansen & Haiman 2003). 

These predictions were strengthened by Miralda-Escude (2003) finding that, 
if the ionizing flux observed today is extrapolated to high z, the photon flux is 
able to reionize the cosmic hydrogen at Zri ~ 8-10 and r ~ 0.08. Furthermore 
Djorgovski et al. (2001) and Becker et al. (2001) reported observations of 
QSO's at 2; > 5, showing an increase of the Lya forest opacity when z goes 
from 5 to 6. Becker et al. (2001) also found Gunn-Peterson effect in one QSO 
at ;2 ~ 6.3. These data conflrmed that a fraction of neutral hydrogen was 
present, in the intergalactic medium, at z ~ 6. (It must be however reminded 
that these very authors warned the reader against an immediate conclusion 
that Zri — 6.) Let us also add that Bullock et al. (2000) revived the idea of the 
excess of galactic satellites, that N-body simulations predict, by considering 
the slowing down of baryon accretion on small velocity lumps after Zri- Their 
analysis apparently leads to requiring Zri < 10, at least for ACDM cosmologies. 

Assuming an (almost) immediate and complete reionization after Zri, in Fig. 1 
we plot how r depends on Zri itself, in various models, assuming Q^h"^ = 0.02. 
A larger r arises in flat models with low Qm, but r > 0.08 can hardly be 
reconciled with Zri < 10 while, for the aforementioned range of Zri, the range 
of T values hes typically below 0.05. 
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Fig. 1. Opacity r as a function of the reionization redshift, assuming instantaneous 
reionization at Zri. Horizontal lines refer to r values 0.08 and 0.16. For r > 0.16 
reionization must have occurred at Zj-i > 15, almost independently of the model. 

Such a satisfactory self-consistent picture has been put in severe danger by 
the recent WMAP release. Using their first-year observation of Ctei at low 
I values, Kogut et al. (2003) estimate that, independently from any model 
assumption, r = 0.17 ± 0.04. However, Spergel et al. (2003) comment that 
the low value of the quadrupole anisotropy and, in general, the rather low 
values of Cti at / < 10 disfavour very high r value, as reionization should 
enhance large angular scales anisotropies (see, however, Efstathiou 2003). In 
fact, once a concordance model is assumed, the likelihood distribution is quite 
fiat towards lower r values and decreases just of 0.05 as r shifts from 0.19 to 
0.11. 

Bearing in mind all these reserves and also the fact that, at the 3-a level, 
T ~ 0.05 is still consistent with observations, if we leave apart any theoretical 
bias, we must also accept that, e.g., a value of r ~ 0.25 is no longer excluded 
by the best observations now available. 

Such findings immediately reopened the theoretical discussion on the reion- 
ization physics. Values of r ~ 0.17 were soon found to be compatible with 
galaxy formation simulations by Ciardi et al. (2003), with reionization caused 
by stars in galaxies with total masses ~ lO^-lO^^M©, assuming a (moderately) 
top-heavy Salpeter IMF, a high stellar production of ionizing photons, and 
a fairly high f^sc — 0.05. These values are however hardly compatible with 
Gunn-Peterson effect findings. 

A more speculative perspective becomes however more realistic in the light of 
high T values: reionization might have occurred twice. Models where this oc- 
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Fig. 2. Dependence of angular power spectra (in K^) on model reionization histories. 
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curred have been proposed by Cen (2003), Wyithe et al. (2003) and Sokasian 
et al. (2003). A rcdshift interval I\z ~ 3, when H is only partially ionized, 
should have elapsed between the two reionization stages. Temporary recom- 
bination might be due to feedback by SNe reducing star formation efficiency 
(see, e.g., Scannapieco et al. 2000, Madau et al. 2001); otherwise, a similar SNe 
action, concentrated in suitable spatial volumes, caused a partial temporary 
recombination in regions where QSO with Gunn-Peterson effect have been 
seen. 



15 
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Fig. 3. Likelihood distributions for the simultaneous determination of reionization 
redshift z^j and ionization rate Xe in models with fixed r = 0.17; the left (right) 
panel corresponds to a 7°-pixel noise cTpix = 2 /xK (0.5 //K). 



Within this context, let us remind that cosmological models with equal r, 
but different reionization histories, predict different CMBR anisotropy and 
polarization. This is shown, first of all, in Fig. 2, where APS, for a fixed r 
value, but different reionization histories, arc shown. The other parameters of 
this spatially flat model are il^/i^ = 0.02, = 0.3, h = 0.7. 

Besides the case of a sharp reionization, taking place at the nearest possible 
redshift, we also consider reionizations starting at larger Zri, followed by a 
single zero-ionization interval, suitably set to recover the fixed r value. Notice, 
in particular, that the high-/, apparently, does not depend on Zri, while the 
shift, at low /, is quite significant. 

The plot in Fig. 2 was made with a suitable generalization of the CMBFAST 
program. Using the techniques described in Sect. 4.1, we can estimate likeli- 
hood distributions and see that the differences shown in Fig. 2 enable us to set 
significant constraints on cosmological parameters, already with a sensitivity 
close to SPOrt's. 

In Fig. 3 (see Colombo et al. 2003; Mainini et al. 2003), we consider two models 
with T = 0.17, either due to full reionization (xg = 1.0) arising at the nearest 
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possible redshift (~ 17) or to a rcionization with constant Xe = 0.6, again 
since the nearest possible redshift (~ 24); the other model parameters are as 
in Fig. 2. 

In both plots 1- and 2-a confidence curves for the simultaneous determination 
of Xe and shown, by solid and dashed lines, respectively. Such curves 

are obtained by averaging over 1000 realizations of each model. The plots are 
also crossed by constant-r curves corresponding to three r values. The l.h.s. 
(r.h.s.) plot corresponds to a pixel noise cXpix = 2 /xK (0.5 /xK), where the pixel 
size is about 7°. 

A first observation is that r is more easily determined than Xe and Zri, sep- 
arately. Furthermore, it is clear that the level of pixel noise that SPOrt is 
likely to achieve (upix — 2 /xK) is just marginal to distinguish between the 
two models. If Upix = 0.5 //K, instead, at the 1-a level, confidence curves are 
disconnected and have just a marginal intersection at the 2-a level. The ca- 
pability of SPOrt to distinguish reionization histories cannot be excluded and 
depends on the model realization in the real Universe. 

In principle, however, the aforementioned figures show that an analysis of 
CMBR polarization data is an excellent probe into the reionization pattern 
of the Universe and, therefore, on birth, evolution and death of early starlike 
objects. 

Let us now come to the point (ii): in shaping polarization spectra, the nature 
of the Dark Energy (DE) component of the Universe plays a role. For the 
sake of example, in Fig. 4 we show the /-dependence of Ctei for / up to 40. 
In all plots Qb = 0.05, h = 0.7, Qm = 0.3, while the solid (dashed) curve 
refers to r = 0.14 (0.20). In all plots the spectra of a ACDM model are shown, 
together with the spectra for models where DE is due to a scalar field 0, self- 
interacting through a Ratra-Peebles (1988, RP hereafter) power-law potential 
V{(f)) = A°'~^^ / . [In some figure of this paper we consider also SUGRA 
potentials \/(0) = (A"+V0") exp(47rG'02)^ g^ax & Martins 1999, 2000]. The 
three plots differ for the choice of the energy scale A and A = logj^ol^/GeV). 
The differences between ACDM and RP models, already significant for A = 2, 
become huge for A = 14, where the first correlation peak has even negative 
values. 

The impact of DE nature can be even better appreciated in Fig. 5, where 
we plot the likelihood distributions on A, as obtainable when artificial data 
built from RP models, with the A scale indicated in top of each frame, are 
"observed". Here again, the techniques described in Sect. 4.1 have been used. 
In the plots, the histograms report the distribution of the peak likelihood, 
as obtained from 1000 realizations of each model. The solid curves, instead, 
are the likelihood distributions, as obtained by averaging over all realization. 
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Upper (lower) plots refer to a pixel noise of 0.2 //K (2.0 //K). 

While, for the unlikely model A = 8, the RP potential is discriminated from 
ACDM at more than S-cr's, at the noise level expected for SPOrt, the discrim- 
ination, with such noise and for lower A's, is not safer than 1 or 2 cr's. These 

plots, however, put in evidence that, even for a = 0.2 fiK, the histogram of 
the peak distribution lays well inside the average likelihood distribution, so 
outlining that we are still far from a regime where cosmic variance dominates. 



A= 14 

1 1 1 1 1 1 1 j 1 1 1 1 1 1 1 1 1 




10 20 30 40 10 20 30 40 10 20 30 40 

1 1 1 

Fig. 4. TE correlation spectra (in K^) for ACDM and RP models compared. Here 
A = log(A/GeV). Solid and dashed lines refer to r = 0.14 and r = 0.20, respectively. 

Besides considering the point (i) and (ii) separately, one should discuss how 
far DE parameters and recombination history can be simultaneously studied, 
at different sensitivity levels. This is an open question that will be debated 
elsewhere. 



All these arguments converge to indicate the importance of inspecting polar- 
ization at low I, where we are still quite far from meeting cosmic variance 
limits. 



2.4 The foregrounds 



Detecting a cosmological signal requires a fair subtraction of all foreground 
signals. Depending on frequency, this may be not so simple; foregrounds, how- 
ever, can be separated from cosmological backgrounds because of their spectral 
properties and their non-Gaussian nature. 

Galactic foregrounds, e.g., have strongly non-Gaussian distributions, so that 
their APS do not provide a complete statistical description. The coordinate 
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X=2 



A=5 



X=8 




Fig. 5. Likelihood distribution on A (solid line), averaged on cosmic and 
noise variance, in measures of artificial data, built from different RP models 
(A = log(A/GeV)). Upper (lower) plots refer to a pixel noise of 0.2 /iK (2.0 fiK). 
Histograms yield the distribution of the likelihood peak in the 1000 realizations 
considered. 
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space description is often more satisfactory. As a matter of fact, the angular 
distribution of total intensity foregrounds has been provided by WMAP. at 
five frequencies between 23 and 94 GHz. The published maps (Bennett et al. 
2003) have a 1° resolution; around this scale, it is recognized that CMBR 
anisotropies dominate in total intensity for Galactic latitudes |6| > 15° and 
frequencies v — 30-100 GHz. 

However, the APS are very convenient for storage of information with a clear 
separation of angular scales, and by analogy with previous CMBR studies, they 
have been applied to foregrounds for about a decade, especially for extensive 
predictions at many scales. Such investigations must be performed at many 
frequencies since the APS of any backgrounds except CMBR depend both on 
/ and V. A common simplifying assumption is 

Cyi = Ay{v)1-'^'' (19) 



with ay independent of both u and /. This is similar to the common real-space 
assumption, that the frequency spectral index be independent of sky position, 
and is likely to be as much unsatisfactory for an accurate description of data. 

Until 2000 the only available experimental knowledge regarded total inten- 
sity, and rather steep APS were often supported for Galactic foregrounds. 
For dust emission a^*^ ~ 3 was found from IRAS (Gautier et al. 1992) and 
DIRBE (Wright 1998); from the combined DIRBE and IRAS maps a slightly 
flatter slope, a^?"* ~ 2.5, was derived by Schlegel et al. (1998). For free-free 
emission, correlating COBE-DMR with DIRBE gave ~ 3 (Kogut et al. 

1996). However modelling free-free emission from Ha maps leads to a flatter 
slope a>p ~ 2.3 and a much lower normalization at 53 GHz (Veeravagham & 
Davies 1997). This discrepancy supported the case for anomalous dust emis- 
sion in the microwave region down to a few GHz, and Draine & Lazarian 
(1998,1999) proposed two possible mechanisms for it (rotational excitation 
of small grains, and magnetic grains). Some further evidence for anomalous 
emission has been proposed over the years [see, e.g. Lazarian & Prunet (2002) 
for a short review] ; however, extensive investigations of available low-latitude 
surveys, carried out within the SPOrt project (Tucci et al. 2000, 2001, 2002), 
supported a more conventional view, that the total intensity APS around a 
few GHz are generally dominated by synchrotron, except in low emission re- 
gions where unresolved point sources dominate especially for / ^ 10^ (cfr. 
Sect. 2.4.2). WMAP's analysis shows that the contribution of spinning dust, 
if any, is less than 5% the foreground intensity at 23 GHz, and this casts 
serious doubts on the relevance of such emission at any frequency. The corre- 
lation of radio and dust emissions is simply explained as due to the fact that 
both trace star-formation activity. WMAP's total foreground APS exhibit 
moderate slopes, ckt — 2 for Z ^ 10^ in all frequency bands. Since thermal 



17 



dust dominates at 94 GHz, we can easily infer a^' ~ 2; the same conclusion 
is likely to apply to free-free and synchrotron which appear to be mixed at 
23-41 GHz. 

For the synchrotron total intensity, results on the slope of the low-frequency 
angular spectrum have been published since 1996, based on several surveys 
such as the 408 MHz maps of Haslam ct al. (1981), the 1420 MHz northern- 
sky survey of Reich & Reich (1986) and the 2326 MHz Rhodes/HartRAO 
survey (Jonas et al. 1998). Slopes ~ 2.5 -j- 3 were given by Tegmark & 
Efstathiou (1996), Bouchet et al. (1996) and Bouchet & Gispert (1999); for 
the Rhodes/HartRAO survey Giardino et al. (2001) report a Galactic latitude 
dependence, with ~ 2.4 in the full maps and ~ 2.9 at |6| > 20° down 
to the resolution limit ~ 1°. However, at the same frequencies an analysis 
of the Tenerife patch reported by Lasenby (1997) provided a nearly scale- 
invariant spectrum aj^ ~ 2; see also the analysis of the 5 GHz Jodrell Bank 
interferometer by Giardino et al. (2000). More recently, local APS slopes in 
the Galactic Plane at 2.4 and 2.7 GHz were found to be ~ 1.7 (Tucci 
et al. 2000; Bruscoli et al. 2002). This result is expected to be affected by 
free-free emitting point sources. However, as already remarked, the WMAP 
full-sky maps give moderate slopes for the APS in all of its frequency bands; 
therefore the low slopes found by Tucci et al. (2000) are confirmed. We can 
infer ~ 2, although a detailed analysis of APS of the individual foreground 
components has not been provided yet. 

Clearly, foregrounds with APS slopes 0:7- ~ 2 or larger are relatively more im- 
portant at large angular scales. An important exception is given by unresolved 
point sources, for which flat intensity APS are predicted, a^?^^ = 0, provided 
their distribution in unstructured (shot noise). Modelling for populations of 
extragalactic radio sources has provided predictions for the amplitude factor 
Ay(z/) (Toffolatti et al. 1998; Guiderdoni 1999); however Galactic Hll regions 
also contribute, and should dominate near the Galactic Plane. Point sources 
must be the most important foreground at sufficiently high /; for instance, for 
the BOOMERANG 153 GHz experiment this happens at the smallest angular 
scales of the measurement, i.e. / ~ 10^ (Netterfield et al. 2002). WMAP's APS 
for / > 200 at 41, 61 and 94 GHz exhibit a similar behaviour, which however 
may be due to instrumental noise. 



2.4.-1 Dust polarization APS 

While waiting for a more complete analysis of the WMAP data, at the time of 
writing direct measurements of foreground Q and U fields, or E and B APS, 
are still lacking. Studies of Galactic polarization APS have been initiated by 
the work of Prunet et al. (1998), who modelled polarized cross-sections as- 
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signing grain compositions and shapes, the distribution of Galactic dust using 
the HI Dwingaloo survey, and the magnetic field by means of gas density 
structures (assuming the field lines to be parallel or perpendicular to gas fil- 
aments). A template for polarized emission from thermal dust was obtained, 
and applying the discrete Fourier transform to it in the small-scale limit, local 
APS were derived for / ^ 2000. The resulting slopes are ue — 1-3, as — 1-4 
and ote — 1-95, and predictions are provided for the scenario to be met by 
PLANCK at 142-217 GHz. Although some assumptions underlying the model 
are now regarded as simphstic by Lazarian & Prunet (2002), the importance 
of this paper should not be underrated, because it clearly recognized that 
polarization APS cannot simply mimic intensity APS up to a constant polar- 
ization degree. In a paper studying the Cui spectrum of starlight (i.e., of the 
distribution of star polarization vectors over the sky), Fosalba et al. (2002) 
notice that polarized emission in the sub-mm/FIR should be related to optical 
polarization by differential absorption, and suggest the use of starlight data 
to improve modelling of polarized dust emission. 

In the light of WMAP and previous results, we do not expect spinning dust to 
be important in polarization. Interesting effects however might arise at some 
frequencies in the magnetic-grain scenario. For grain models with random Fe 
inclusions the polarized emission appears to be peaked around v ^ 80 GHz, 
and differing from the usual expectation for thermal emission (but also for 
spinning dust), the polarization vector should be parallel to the Galactic mag- 
netic field, as is the case for starhght polarization by differential absorption. 
APS modelling has not been provided so far for the anomalous dust emission; 
the best estimates for spinning dust provided by Tegmark et al. (2000) simply 
assume the same APS behaviour as for thermal dust. 



2.4-2 Galactic synchrotron APS 

Differing from total emission, polarized Galactic emission is not expected to be 
much affected by free-free in the microwave region. Most investigations have 
thereby been focused on synchrotron. Starting with the work of Tucci et al. 
(2000), a number of papers tackled the problem of polarized synchrotron APS 
at low frequencies in view of extrapolations to the cosmological window. In 
particular, within the ambit of the SPOrt project the following surveys were 
analysed: 

• The Parkes 2.417 GHz survey of the southern Galactic Plane, described by 
Duncan et al. (1997) (henceforth, D97). It covers a strip, 127° long, with a 
FWHM resolution of 10.4'. 

• The Effelsberg 2.695 GHz survey (Duncan et al. 1999: D99), covering a 70° 
strip of the first-quadrant Galactic Plane with a resolution of 5.1'. 

• The Effelsberg 1.4 survey (Uyaniker et al. 1999: U99) with five disjoint 
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regions fairly close to the Galactic Plane, —15° < b < 20°. The angular 
resolution is there 9.35'. 

• The Leiden survey of Brouw & Spoelstra (1976) (BS76). It covers about 40% 
of the sky up to near the north Galactic Pole at five frequencies, between 
408 and 1411 MHz, with resolution varying between 2.3° and 0.6°, but with 
a significant undersampling (BS76). 

• The ATCA intcrfcrometric survey of the SGPS Test Region, a 4° x 7° region 
internal to the Parkes survey, with an elliptic beam of 67" x 87" (Gaensler 
et al. 2001: GOl). Only the 1.404 GHz channel has been used so far for 
computation of APS. 

The D97 and D99 surveys showed that Galactic-Plane polarization is not much 
correlated to total intensity, being often high in the absence of sources in total 
intensity, and exhibits a foreground component generally supposed to originate 
locally, i.e. within a few kpc. Local APS have been computed in twelve D97 
patches and six D99 patches of size 10° x 10° (Tucci et al. 2000, 2001) and ex- 
hibit moderate slopes — as — 1-5, with no systematic differences between 
high- and low-polarized emission regions. A recent re-analysis of D97 data by 
another group (Giardino et al. 2002) gives similar results. Out of the Galactic 
plane, the U99 survey shows larger dispersions in the spectral parameters, due 
to structures with strong polarized signals (e.g., in the Cygnus and Fan re- 
gions) and regions of very low total intensity. Five low-intensity patches with 
ax — ^1 which must be dominated in intensity by discrete sources, have been 
used to infer a rather accurate estimate C|^^'' ~ 5 x 10~^ for the sum of 
all point source populations at 1.4 GHz. Adopting a radio-source polarization 
degree ~ 5%, the estimate C%^^ ~ 1.3 x 10"^ was provided at 1.4 GHz, as 
well as other estimates at different frequencies with reasonable assumptions on 
the frequency spectral index (Bruscoli et al. 2002). Such estimates imply that 
the contribution of point sources should be negligible for all of the polariza- 
tion APS reported for all of the aforementioned surveys, as confirmed by the 
fact, tliat the slopes of polarization APS do not show any fiattening towards 
higher /. The impact of Faraday rotation, however, is not clear from the data, 
and must be investigated looking at different frequencies and angular scales. 
A multiwavelength study has been performed in the range / < 70 by Bruscoli 
et al. (2002) using three large patches of the BS76 survey at various latitudes 
up to near the Galactic Pole. At such scales no definite trends were found 
in the spectral slopes for increasing Galactic latitude. A clear fiattening of 
polarization APS, aE,B — 1, appears at low frequencies, u < 610 MHz, where 
Faraday rotation is larger; on the other hand, at 1.4 GHz no significant depen- 
dence on Galactic latitude appears, due to the large error bars. Also the APS 
normalizations are essentially consistent, as shown by Fig. 6. According to the 
analysis of D97, D99 and BS76 data, the Cm spectra at z/ > 1.4 GHz may be 
only slightly steeper than Cpi, since the best slope estimate is an — 1.6 -r- 1.7 
and the difference has little statistical significance. The substantial agreement 
does not prove, however, that the polarization APS at > 1.4 GHz do reflect 
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Fig. 6. The APS for E mode (thick and dash-dotted hnes) and polarized intensity 
(thin hnes). The CMBR concordance model assumes the parameters adopted by 
Kovac et al. (2002). The thin and dash-dotted CMB lines refer to a model with 
a reheating optical depth r = 0.1; the polarized intensity curve is a power-law 
approximation to the output of simulations reported by Tucci et al. (2002). The 
synchrotron APS arc best-fit curves, extrapolated to 30 and 90 GHz from the 
following maps: The template of Bernardi et al. (2003), / = 2-20; olation; the 
intermediate-latitude patch in BS76 (Bruscoli et al. 2002), I = 10-70; the D97 
survey (Tucci ct al. 2000), / = 100-800. For the template of Bernardi et al. (2003) we 
also give the APS from reduced maps at high Galactic latitudes, with extrapolations 
to / = 20-1000 represented by dotted lines. 
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the intrinsic properties of synchrotron at all scales. 

A more recent analysis of Tucci et al. (2002) in fact shows that this is not the 
case at 1411 MHz at high resolutions. Rich small-scale structures appear in 
the ATCA maps (GOl) and other polarimetric observations. The sharp-edge 
patches in the GOl maps of polarization angles are particularly impressive, and 
strongly suggest that many structures are in fact produced by sharp Faraday- 
rotation changes in the polarization angle. This has a strong bearing on Cei 
and Cbi'- The polarization APS computed in a 4° x 4° box and several smaller 
patches give aE — cub — 2.8, whereas an — 1.6. Such values have been com- 
pared to those computed in an overlapping 5° x 5° box extracted from the D97 
survey. Large discrepancies are found for aE and a^, Aa ~ 1, but not for an- 
Rescaling the data with a synchrotron spectral index /Jgync between —2.5 and 
—3, The ATCA Cnz spectrum appears to be the natural extension to smaller 
scales of the corresponding spectrum at I < 800, since the normahzation is also 
quite consistent. On the other hand, the power excess of E and B modes in 
the ATCA spectra is nearly one order of magnitude at scales 600 ^ / ^ 1000 
and declines at / greater than a few thousand. This excess must be due to 
Faraday boosts of the polarization angles (Tucci et al. 2002). 

2.4-3 A synchrotron template at SPOrt scales 

In the light of the aforementioned results, a simple extrapolation of Cei 
and Cbi from ~ 1 GHz to higher frequencies is not possible using an l- 
independent spectral index. How can we get rid of Faraday rotation? Cia- 
rdino et al. (2002) use D97, the total-intensity survey of Haslam et al. (1981) 
and other radio data to produce a real-space model of polarized synchrotron 
emission T and H; the polarization angle however is just taken as a random 
quantity, so the final synthetic maps were declared to be a toy model. A more 
powerful approach has been recently introduced by Bernardi et al. (2003), 
intended to produce a Faraday-free template of polarized synchrotron. It is 
based on the following steps: (a) A synchrotron intensity map is built using 
the 408 MHz total intensity survey of Haslam et al. (1981) and the 1.4 GHz 
survey of Reich & Reich (1986). Two different frequencies are necessary to 
clean the maps from free-free emission, (b) The polarized intensity H is mod- 
elled, assuming that polarization is produced within a spherical "polarization 
horizon" and calibrating a normalization factor with comparison with BS76 in 
the Fan region, (c) A map of polarization angles is built from starlight polar- 
ization data, with interpolation from Heiles (2000) catalog and a 90° rotation. 
The last step, which assumes an essentially common location for dust selective 
absorption of starlight and for synchrotron emission, is of crucial importance. 
Several tests discussed by Bernardi et al. (2003) support the reliability of the 
template under this respect, and in connection with the polarization horizon 
modelling for H as well. Due to the Heiles data sampling the template has 
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been provided for an angular scale of 7°. An update of the template, taking 
advantage of the WMAP intensity maps at 23 GHz, is in preparation. APS 
computed on Faraday-free templates can be extrapolated to higher frequen- 
cies better than those of current surveys. Due to the low resolution, results 
could be provided by Bernardi et al. (2003) only for I < 20, where the full-map 
template gives ap = 1.9 ± 0.2 and an = 1-9 ± 0.2. The substantial agreement 
with surveys at 2.4-2.7 GHz reinforces both the validity of the template and 
our previous support to moderate spectral slopes. The high slopes aE,B — 2.8 
found by ATCA in the SGPS Test Region at 1.4 GHz, on the other hand, must 
overrate intrinsic synchrotron slopes substantially, being due to the high-Z de- 
cline of the spectral power excess. 

In Fig. 6 we have collected the extrapolations to 30 and 90 GHz of some of the 
power-law fits to £'-mode and polarized intensity APS, without attempting to 
give error bars. Such extrapolations, which at present provide only reasonable 
guidelines, cover the entire range / = 2-1000 including results from various 
surveys at z/ = 1.4 -h 2.4 GHz and the aforementioned template. For the tem- 
plate we show the APS derived from both full maps and reduced maps (50%) 
containing only the lowest-signal pixels; for the latter, the extrapolation to the 
angular range I — 20-=- 1000 is also provided. Clearly, in view of the differences 
in the investigated sky regions, coverages and frequencies, one should not ex- 
pect to find a detailed agreement. However, the curves displayed in Fig. 6 are 
roughly consistent with each other even for the normalization factor. Only 
the template APS derived from the reduced maps (with low-signal pixels) 
have a lower normalization, and suggest that the detection of the secondary- 
ionization spectral bump should not be hampered by synchrotron at 90 GHz. 



3 The SPOrt Experiment 



The SPOrt experiment, carried on under the scientific responsibility of an 
international collaboration of Institutes headed by the lASF-CNR in Bologna, 
and fully funded by the Italian Space Agency (ASI), was selected by ESA in 
1997 to be fiown on board the International Space Station (ISS) during the 
Early Utihzation Phase. 

The payload, shown in Fig. 7, houses four corrugated feed horns feeding a 
22, a 32, and a pair of 90 GHz channels providing direct measurements of 
the Q and U Stokes parameters. It will be located on the Columbus External 
Payload Facility of the ISS in 2006, for a mission with a minimum lifetime of 
18 months. 

The SPOrt antennae, looking at the zenit, will cover 80% of the sky thanks to 
the motion of the Space Station, whose orbit is tilted by 51.6° with respect to 
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Table 1 

SPOrt main features: A^pix is the number of FWHM pixels covered by SPOrt, ais 
is the istantaneous sensitivity (1 second), and cjpix is the pixel sensitivity for a 
two-year mission. 

I' channels BW FWHM Orbit Time Coverage A^pix o"is Cpix 

(GHz) (#) n (s) (%) (mKsV^) (^K) 

22 1 0.5 1.6 

32 1 10% 7 5400 80 660 0.5 1.6 

90 2 0.57 1.8 
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Fig. 8. Left: The sky, in Celestial coordinates, as scanned by SPOrt in few orbits. 
Right: Pixel observing time (seconds) for two years of data taking. The pixel size is 
about 7° (HEALPix parameter Nside = 8). 

Table 2 

SPOrt sensitivity to CMBR, in thermodynamic temperature, combining the four 
channels and including Galactic synchrotron subtraction, for a two-year mission. 





"^pix 


'^(Prms) 


(mKsi/2) 


(mK) 


(^K) 


0.53 


1.7 


0.15 



the Celestial equator and is characterized by a period of about 5400 s and a 
precession of about 70 days. This results in the sky scanning pattern shown in 
Fig. 8, providing a full coverage of the accessible region every 70 days. In Fig. 8 
a map of the time spent over each pixel, of about 7° (HEALPix^ parameter 
Nside = 8), is also shown. 

The main features of the SPOrt experiment, summarised in Table 1, have 
been chosen to allow the accomplishment of SPOrt primary goals, e.g. a ten- 
tative detection of CMBR polarization on large angular scales, and the map- 
ping of Galactic synchrotron emission. The former asks for both a nearly all- 
sky survey and a multifrequency approach to control possible contaminations 
from Galactic foregrounds; the latter is expected from the channels at 22 and 
32 GHz. 

The direct and simultaneous measurement of both Stokes parameters Q and U 
optimizes the observing time efficiency, a remarkable improvement compared 
to other schemes providing either Q ot U and thus reducing it by a factor 2. 

The SPOrt expected sensitivity to CMBR, resulting from the combination of 
the four channels and taking into account the worsening due to foreground 
subtraction (Dodelson 1997), is reported in Table 2. In the frame of recent 
WMAP results on the optical depth of the Universe at the re-ionization epoch, 
SPOrt promises a solid detection of large scale CMBR polarization. 

^ http:/ /www.eso.org/science/healpix 
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Fig. 9. Schematic block diagram of the SPOrt radiometers. Polarizer and OMT ex- 
tract the two circularly polarized components LHP and RHP collected by the horn. 
After amplification, the correlation unit (based on a Hybrid Phase Discriminator, 
see Peverini et al. 2002 for details) provides directly both Q and U. 



3. 1 Instrument Design and Analysis 



The main concern in designing an experiment to measure CMBR polarization 
is the low level of the expected signal (1-10% of the already tiny temperature 
anisotropies, depending on the scale), requiring specific instrumentation. 

The design optimization with respect to systematics generation, long term sta- 
bility and observing time efficiency is indeed even more critical here than for 
experiments primarily designed to investigate CMBR temperature anisotropies. 
This is demonstrated by previous attempts to measure CMBR polarization 
where the final sensitivity was sizeably worse than that expected from the 
istantaneous sensitivity of the front-end noise (Keating et al. 2001; Hedman 
et al. 2002; Kogut et al. 2003). 

The following main choices were made while designing the polarimeter: 

• correlation architecture to improve the stability; 

• correlation of the two circularly polarized components and to directly 
and simultaneously measure both Q and U from the real and imaginary parts 
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of the correlation, respectively, 

Q oc ^{ErED ; U oc "^{ErED (20) 

while keeping 100% efficiency; 

• detailed analysis of the correlation scheme to identify critical components 
and set their specifications in order to keep the systematics at a level suitable 
to measure CMBR polarization; 

• custom development of components when state-of-the-art is not compliant 
with SPOrt requirements; 

• on axis simple optics (corrugated feed horns) to minimize the spurious po- 
larization induced by both the / pattern (see Carretti et al. 2001 for its 
definition) and the CMBR temperature anisotropy on the beam scale. 

The equation giving the radiometer sensitivity (WoUack 1995; WoUack et al. 
1998) 

/t.2 / /\r\ 2 

= + T„2ff.et (^) + A7?ff.et (21) 

can help us find the parameters to be controlled to minimize the systemat- 
ics. Here Tgyg, Toffset and AToffset represent the system temperature, the offset 
equivalent temperature and its fiuctuation, respectively; G is the radiometer 
gain, t the integration time, Az/ the radiofrequency bandwidth and A; ~ 1 a 
constant depending on the radiometer type. 

The first term in Eq. (21) represents the white noise of an ideal and stable 
radiometer. Its contribution can be reduced by cooling to cryogenic tempera- 
ture the front-end (first stage of the Low Noise Amplifiers (LNAs), polarizer 
and the OMT). Unfortunately, passive cooling is not feasible for low orbits 
because of variations of the Sun illumination. For SPOrt, a mechanical pulse- 
tube cryocooler was adopted, ensuring high cooling efficiency down to ~80- 
90 K. The needed temperature stability (±0.1 K) is achieved by a closed-loop 
active control. An active control is also adopted for the warm parts (Horn and 
back-end), allowing a temperature stability of ±0.2 K. 

The second and the third term of Eq. (21) represent the contributions of gain 
and offset fiuctuations, respectively, generated by instrument instabilities. It 
is clear that the noise behaviour can be kept close to the ideal (white) case 
provided the offset is kept under due control. 

A scheme of the SPOrt radiometers is sketched in Fig. 9. The Polarizer and the 
OMT extract the two circularly polarized components El and Er collected 
by the dual-polarization feed horn. After amplification, the two components 
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are correlated by the Correlation Unit (CU), consisting of an Hybrid Phase 
Discriminator (HPD), four diodes and two differential amplifiers. In details, 
the HPD generates the four outputs: 



Ci(xEl + Er (22) 

C2^El-Er (23) 

Cs^El+jEr (24) 

C4(xEl-jEr (25) 

which are square-law detected by the four diodes 

\C,f oc IElI' + \Er\^ + 2^{ErEI) (26) 

IC2I' oc \ELf - IEr^ - 2^{ErEI) (27) 

IC3I' oc \El\^ + \Er\^ + 2%ErEI) (28) 

IC4I' oc \El\^ - \Er\' - 2'^{ErEI). (29) 



The differences performed by the two differential amplifiers provide 

\C,\^ - \C2f oc ^{ErED Q 
IC3I' - IC4I' oc Q{ErEI) oc U 



allowing the simultaneous measurement of the two Stokes parameters Q &i U. 
A nice feature is that no effort is needed to equalize the mean phase difference 
between Er and El when propagating the signal through the instrument, since 
an error on this difference simply results in a polarization angle rotation, well 
recoverable in the calibration procedure. Indeed, linear polarization might also 
be measured by correlating the two linear polarizations E^, Ey. However, this 
would only provide one linear Stokes parameter (U), reducing the experiment 
efficiency, and would generate a combination of U and V in the output in 
case of errors in the phase difference. Therefore, great care would be needed 
in equalizing the path of the two components. 

Our analysis to track down critical components for the offset generation pointed 
out important sources in both the CU and the antenna system (horn, polarizer 
and OMT). 

We found that the CU needs an HPD with high rejection of the unpolarized 

component, at a level well beyond available instrumentation. A custom device 
(Tascone et al. 2002; Peverini et al. 2002) providing > 30 dB rejection was 
thus developed. In combination with a lock-in system, this makes the CU 
contribution to the offset negligible, the total rejection being ~ 60 dB. 
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The contribution to the offset coming from the antenna system is given by 
(Carretti et ah 2001): 

(rppol\ 
Tsky + T^°[^ ^ 1 , (31) 



where Tsky is the signal coming from the sky, T^^°{^ is the noise generated by 
the horn alone, T^^ise noise temperature of the whole antenna system, r] 

is the efficiency of the feed horn and T^^^ is the physical temperature of the 
polarizer. The two quantities 
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describe the performances of the OMT and the polarizer, respectively, in terms 
of offset generation. Uncorrelated signals like noise and sky emission are par- 
tially detected as correlated because of the OMT cross-talk {Sai and Sbi 
representing the transmission and cross-talk coefficients of the OMT, respec- 
tively) and the polarizer attenuation difference {S\\ and S± giving the trans- 
missions of the two polarization states). 

To set specification requirements for our OMT and polarizer we considered the 
radiometer knee frequency, /knee- This parameter provides the time scale at 
which the low frequency component of the noise, also known as 1//, prevails 
on the white one, and is therefore suitable to quantify radiometer instabilities. 

As is known, destriping techniques can remove most of the effects of 1// noise 
provided the radiometer knee frequency is lower than the signal modulation 
frequency (see Sect. 3.3) that, for SPOrt, corresponds to the orbit frequency 
/o = 0.18 mHz. 

Currently available InP LNA have rather high knee frequencies (/k^ee ~ ^00- 
1000 Hz). However, the knee frequency of a correlation receiver is related to 
that of its amplifiers by the formula: 

Anee=f^)'''"/,^n:e (34) 



where ToSset is the radiometric offset, Tgys the system temperature and a ~ 1 
(WoUack 1995) is the power-law index of the LNA's 1// noise behaviour. 



29 




Fig. 10. Magnitude of both the isolation between the two rectangular ports and the 
cross-talk between the two polarization states for the 32 GHz OMT of SPOrt. The 
vertical dotted lines define the 10% band. 
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Fig. 11. Co-polar, cross-polar and / patterns, with respect to the axis distance 9, 
for the 90 GHz feed horns of SPOrt . The spurious polarization coming from the 
horn is generated by the anisotropy distribution l\T{0(j)) in combination with the 
/ pattern (see text). 
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To match the condition for a successful destriping (/knee < /o, with /knee = 
O.l/o = 0.018 niHz as a goal), SPOrt requirements on the OMT cross-talk 
and the difference between the attenuations of the two polarization states were 
quantified in —60 dB and —30 dB, respectively. This guarantees an offset value 
as low as ToSset ~ 50 mK which, combined with a Tgys ~ 100 K, translates 
into the knee frequency 

/knee - 2.5 X lO'^/^^fee (35) 



fully satisfying the requirements for an efficient destriping. 

Commercially available statc-of-thc-art OMTs have isolation values not bet- 
ter than —40 dB, which forced us to develop new hardware. Fig. 10 shows the 
extremely good results already obtained for the 32 GHz channel: the cross-talk 
is as low as —65 dB, with an isolation of about 70 dB. 

Besides the offset generation, a further worrysome source of systematic errors 
is the spurious polarization generated by the optics (Carretti et al. 2001). This 
is due to the anisotropy distribution of the unpolarized radiation, modulated 
by the / pattern: 



AT 



ihorn 



1 



TT n/2 





- AT,(^, + 7r/2) + AT,{e, + tt) 

- AT6(^,0 + 3/27r)] f{,eA)dct>, 



(36) 



+ xi0,cf>)P*i0,<P + 7r/2), (37) 

where P and x are the co-polar and cross-polar patterns, normalized to the P 
maximum, respectively, and Qa is the antenna beam. The / pattern consists 
of a 4-lobe structure, with lobe size close to the FWHM of the instrument. A 
worst-case analysis of the contamination gives 

Ar^°- = SPhorn ATr„,,(FWHM) , (38) 



with 



7r/2 

SPhorn-2 — |sin^d^ J d<Pf{9,<P) 





(39) 



and Arrms(FWHM) the rms temperature anisotropy on FWHM scale. As 
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shown in Fig. 11, in case of the SPOrt feed horns the contribution of the 
/ pattern is 5Phorn ~ ~24 dB and the rms contamination due to the 30 /xK- 
CMBR anisotropy is thus lower than 0.2 /iK. 



3.2 Calibration procedure 



The accuracy needed for measuring CMBR polarization requires good meth- 
ods for calibrating the response of the instrument to small polarized signals. 
Furthermore, in the absence of well-characterized astrophysical sources, spe- 
cialized techniques are needed to inject calibration markers. In fact, standard 
marker injectors are not suitable for calibrating tensorial quantities as the 
{Q, U) pair measured by SPOrt. A new concept calibrator, valid for any radio- 
polarimeter and based on the insertion of three signals at different position 
angles, has thus been developed (Baralis et al. 2002). 



To focus the problem it is convenient to represent the entire radiometer in 
terms of a system with two input and output signals. The two inputs A and 
B are high frequency signals in the circular polarization base {El and Er)^ 
whereas the two outputs and Um are the low frequency signals correspond- 
ing to the Stokes parameters directly measured by the radiometer. According 
to this model, the output of the radiometer can be described by the following 
matrix expression: 
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(40) 



where Q and U are the input Stokes parameters. The matrices ^ and ^ 
in Eq. (40) are generic, but real, since they deal with quadratic quantities. 
Moreover, they transform polarization circles in the QU-plane into rotated 
and translated ellipses. Under the reasonable assumption that the intensities 
of the total signals A and B are equal (|Ap = \B\'^ oc P), Eq. (40) can be 
rewritten as Yi — ^2L + QP, where C is a two-element column vector. The 
calibration procedure, i.e. the evaluation of the matrices H_ and the vector 
C_, can be accomplished by directly measuring the quantities Qm and Um 
in presence of predefined signals (markers). In fact, by injecting at different 
times three markers with known polarization angles 9i, each corresponding to a 
linearly polarized field, and by detecting the corresponding output variations, 
one ends up with a matrix equation to be solved for the matrix 

H-rMa _ M1M2 _ Ml 1 = _ Ml M.Mii an 

= 1- AP3 APi AP2 APi J I- AP3 APi ' AP2 APi J V ^ 

where AXj, AFj and APj are the variations of the input and output Stokes 
parameters and of the total power, respectively. To obtain a well-conditioned 
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coaxial port 




(b) 



Fig. 12. Scheme of the marker injector, a) Longitudinal cut. b) Cross-section of the 
central block. 

matrix the markers should correspond to signals with a relative rotation of 
either 45° or 60°. 



3.2.1 Marker injector 

The calibration procedure described in the previous section requires the injec- 
tion of different markers in the radiometer, just behind the antenna. Unfor- 
tunately, this procedure must be carried out in operative conditions, e.g. in 
presence of both polarized and unpolarized radiation. Therefore, the marker 
injector must degrade the relevant signals as little as possible. In particular, 
this device has to exhibit a high return loss and prevent depolarization of the 
incoming signals. 

The design of the Ka-band prototype is depicted in Fig. 12, and consists of 
three blocks (Peverini et al. 2002b). The central block is formed by an enlarged 
circular waveguide in which eight K-connectors are inserted. Although only 
three markers have to be injected inside the circular waveguide, all the eight 
connectors are necessary to preserve the azimuthal symmetry of the structure 
and, hence, to avoid depolarization of the incoming signals. The coupling level 
between coaxial and circular ports can be controlled by adjusting the pene- 
tration length of the internal wire of the K-connectors. The two lateral blocks 
provide a matching structure between the input/ouput circular waveguides 
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Fig. 13. Measured an computed scattering parameters of the Ka-band marker injec- 
tor. Top: Reflection coefficient. Bottom: Coaxial/circular coupling factor. Vertical 
dashed lines define the SPOrt band. 
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and the enlarged central waveguide. The matching section is formed by a cas- 
cade of several circular steps. The device was analyzed by the Generalized 
Scattering Matrix approach (Itoh 1989) and designed by evolution-strategy 
methods. All the components were manufactured by highly precise electrical- 
discharge techniques. 

In Fig. 13 the measured and simulated reflection coefficient 5*11 at the circular 
port for the fundamental TEn mode and the coupUng 5*31 between circular 
and coaxial ports are reported. The measured coaxial/circular coupling factor 
is ~ —42 dB in all the band of interest (30.4 ^ 33.6 GHz) and the fundamental 
mode at the circular port exhibits a reflection coefficient better than — 35dB 
in the same band. The measured in-band insertion loss is about 0.02 dB, 
which can be correctly predicted by using an equivalent surface resistance of 
10 nflcm. 

3.3 Destriping and Map Making 

Low frequency noise induces correlations among successive samples of the 
measured signal and leads to typical striping effects when producing sky maps. 

However, when data are taken from spinning spacecrafts, most of the low 
frequency noise can be removed by software, provided the radiometer knee 
frequency is lower than the satellite spin frequency (Janssen et al. 1996). As 
detailed in Sect. 2.1.1, the low frequency noise of the SPOrt radiometers is ex- 
pected to have a l//-like spectrum dominated by transistor gain fluctuations. 
The present state of the instrument already guarantees a knee frequency /knee 
lower than the ISS orbit frequency /q. 

Various destriping algorithms have been proposed in recent years to clean 
CMBR anisotropy data (e.g. Delabrouille 1998; Maino et al. 1999), together 
with a first algorithm specifically designed for the polarization case (Revenu 
et al. 2000). They are generally based on minimizations and involve large 
matrix inversions. 

A different technique has been implemented for SPOrt (Sbarra et al. 2003), 
consisting of a simple but effective iterative algorithm relying upon minimal 
assumptions: the radiometer must be stable during the signal modulation pe- 
riod (the time needed to complete one orbit in case of SPOrt), so that the 
instrumental offset can be assumed to be constant over the same period, and 
there must be enough overlap between different orbits. In such a situation the 
noise can be split into two parts: for time scales shorter than the orbit period 
it is essentially white, whereas for longer timescales the 1/f component can 
be modelled as a different constant offset for each orbit. A simple iterative 
procedure can then be applied to remove these offsets from the Time Ordered 
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Fig. 14. Noise simulated maps before (top) and after (bottom) destriping, for the 
case /knee = 1-8 • 10^, in Kelvin. The HEALPix parameter Nside is 32 (courtesy of 
A&A). 

Data (TOD) before map-making. No assumptions need to be made about the 
statistical properties of the noise. The map-making itself consists of a simple 
average of the measurements corresponding to the same pixel, which is the 
optimal method once only white noise is left (Tegmark 1997). 

Although the algorithm has been studied to destripe Q and U data, it can be 
easily simplified to deal with scalar quantities. The average of the measured 
signal is lost in the latter case. However, for maps of polarization data hke Q 
and [/, the average signal can still be kept provided the polarimeter reference 
frame rotates about the standard reference frame (polar basis, see Berkhuijsen 
1975) while nmning along each orbit, as in case of SPOrt. This is a nice feature, 
especially when measuring foreground contributions. 

The good performances of the technique are made evident in Fig. 14 where we 
show a simulated noise map (both 1// and white) before and after destriping. 
One possible way to quantify the quality of the destriping is measuring the 

Table 3 

Excess rms noise due to low frequency contributions, with respect to the white noise 
level, for pixels of ~ 7" . 



/knee (Hz) 


Before Destriping 


After Destriping 


1.8 X 10-^ 


310% 


6% 


1.8 X 10-5 


35% 


< 1% 



fractional excess pixel noise with respect to the case of purely white noise. 
Results are shown in Table 3 for two different values of the knee frequency, 
corresponding to the SPOrt goal knee frequency, /knee = 1-8 x 10"^ Hz, and 
the SPOrt orbit frequency, /o = 1.8 x 10~^ Hz, the latter representing a 
conservative case. 

Another way to quantify the residual correlated noise is measuring and in- 
specting the two-point correlation functions C^{0) = {Q{1)Q{2)) and C'" {6) — 
(JJ{1)U{2)) (see Sect. 2.1.1) of simulated Q and U noise maps. Averages and 
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Fig. 15. Average and la band of correlation functions C^{9) measured from 500 
noise maps, before and after destriping, for two values of the knee frequency /knee- 
Before destriping the y scale is not the same for the two /knee values. Each simulation 
corresponds to about one year of SPOrt realistic data taking. The case of purely 
white noise is shown for comparison (courtesy of A&A) . 
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Fig. 16. Increment, due to the presence of residual correlated noise, in the rms 
measured noise power spectra, after destriping, as a function of the multipole, 
percentage of the rms of purely white-noise power spectra (courtesy of A& A) . 
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Fig. 17. Synchrotron template at 22 GHz (see Sect. 2.4.3), featuring a Prms = IT^K. 

la bands of 500 correlation functions C'^(^) measured from maps containing 
both white and 1// noise, before and after destriping, are compared to the 
purely white noise case in Fig. 15 for the same knee frequencies as in the 
previous test. As expected, the correlated noise is strongly reduced by the 
destriping procedure, the residuals falling within the statistical error of the 
white noise case for /knee = /knee- 

The polarization power spectra Cei and Cbi can be obtained from the mea- 
sured correlation functions by inverting Eqs. (9) (Kamionkowski et al. 1997b; 
Zaldarriaga 1998; Ng & Liu 1999). If the correlation functions are measured 
directly on Q and U maps via 0{Np^^) operations (iVpix being the number of 
pixels in the measured map), this method (Sbarra et al. 2003) has the ad- 
vantage of avoiding edge problems (Chon et al. 2003) arising when using fast 
spherical transforms. 

The region of low multipoles is the most sensitive to low frequency residuals, 
some contributions being always found here even after the application of other 
destriping techniques (Maino et al. 1999). Even though residual noise corre- 
lations can be modelled and subtracted from the measured C^'^{9) functions 
before integration, their presence translates into an increment of the rms of 
the measured power spectra. A rough estimate for the case of SPOrt is shown 
in Fig. 16. 



4 Astrophysics and Cosmology with SPOrt 

A first goal that SPOrt will achieve is a deeper knowledge of the polarization 
of the Galactic signal, in the microwave band. In particular, according to the 
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Fig. 18. 1 — a error bands, as expected from SPOrt measurement of the Cei power 
spectrum, including cosmic variance, considering the sensitivity quoted in Table 2 
and r = 0.17, as recently found by WMAP. APS corresponding to r = 0.1 and 
T = 0.05 are shown for comparison. 



synchrotron template described in Sect. 2.4.3, which we believe to be not 
too far away from reality, SPOrt should be able to produce polarized Galaxy 
maps at both 22 and 30 GHz. In fact, as shown in Fig. 17, the synchrotron 
polarized emission, as seen by the 7°-FWHM SPOrt optics at 22 GHz, is 
expected to have a peak of about 100 /iK, two orders of magnitude larger 
than the experimental pixel sensitivity in the same channel (o"pix = 1-6 /iK); 
furthermore, the average of the polarized signal is Prms = 17 /iK. Assuming a 
power-low dependence on the frequency with index /3synch = —3, the polarized 
intensity at 30 GHz is then expected to be just a factor 2-3 lower than at 
22 GHz, so that SPOrt is likely to provide a Galaxy map at this frequency as 
well. 

The SPOrt pixel sensitivity quoted in Table 2 does not allow the building of 
CMBR polarization maps. However, a CMBR band power spectrum as that 
shown in Fig. 18 is expected to be provided if the recent WMAP results are 
confirmed, l-o" error bands are calculated following the recipe provided by 
Tegmark (1997). The plot also suggests that a non-null (at 95%C.L.) mea- 
surement of the mean polarized signal Prms might be provided, already for 
T ~ 0.1, by full-sky statistical analyses like the Maximum Likehhood flat- 
spectrum analysis developed by Zaldarriaga (1998). 

In addition to measuring two low-multipole bands of the power spectrum, 
combining SPOrt polarization data with WMAP temperature data will pro- 
vide a measurement of the Cf^ power spectrum where errors on T will be 
completely decorrelated from those on Q,U. 
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4-1 Measuring the reionization history of the Universe 

In this subsection we briefly describe our procedure to work out cosmological 
parameters from SPOrt data, stressing, in particular, our expected capacity 
to inspect the actual value of the cosmic opacity r to CMBR photons. Q and 
U maps are assumed to be free from foreground contamination, the sensi- 
tivity worsening coming from this separation, estimated with the procedure 
suggested by Dodelson (1997), being already included in the quoted pixel sen- 
sitivity (see Table 2 in Sect. 3). 

In particular, we shall recall and extend previous results by Colombo & Bonomctto 
(2003), showing how measurements on r and Us (primeval spectral index) are 
far less degenerate when polarization data are exploited, at variance from what 
happens if anisotropy data only are used. 

SPOrt large beamwidth allows us to bin data in a fairly low number of pix- 
els (~ 2400 for a HEALPix parameter Nside = 16, yielding a mean an- 
gular distance ~ 3.5° between pixel centres); therefore, working in the an- 
gular (pixel) space is not numerically expensive; moreover this choice eases 
the combination of data from experiments with different sky coverages and 
detector speciflcations, as we expect to have for anisotropy and polariza- 
tion: anisotropy data specify T on Nt pixels and polarization data specify 
Q and U on Np pixels. We order data into a {Nt + 2A^p)-component vector 
X = {T{i = 1, Nt), Q{i = 1, Np),U{i = 1, Np)}. 

We build artificial data sets, binned according to HEALPix, by using CMB- 
FAST. Data are realizations of a given cosmological model Ai and, in each 
data-map, a white noise map in included. Therefore, the correlation matrix 
(x'^iXj) = Cij — Sij + Njj (brackets mean ensemble average) is the sum of a 
signal term S^j, which depends on cosmology and can be evaluated using the 
relations given in Sec. 2.1.1, and a term Njj, which accounts for the detector 
noise. We assume no correlation among noise in different modes and pixels, so 
that Njj is diagonal and has distinct values erf, and Up, in Nt and 2A^p pixels, 
respectively. 

Making use of the matrix Cjj, the likelihood of a model A4' given the synthetic 
data X reads: 

£(A^'|x) = [(27r)^^+'^^ det C']-^/' exp[-x'^C'~'x/2] . (42) 

Confidence regions in parameter space are found by a likelihood ratio criterion. 
This procedure is applied to extract cosmological parameters, by inspecting 
sets of models M.' against real data. 
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Confidence regions in the n^-r plane are obtained while keeping the other 
parameters at the fixed values: = 0.35, Q\ = 0.65, = 0.05, h = 0.65, 
and assuming no tensor modes. Disentangling r and measurements is a 
typical degeneracy removal allowed by polarization data. 

All models considered here reionize at the nearest possible redshift, where 
baryonic materials undergo a (nearly) sharp transition from — to Xe — 1- 
In Sec. 2.3 we outlined that this simple reionization history, not only is not 
unique, but could also be rather unlikely, if the highest values of r found by 
the first-year WMAP release are confirmed. In that section, however, we also 
showed that detecting details in the reionization history, at our expected level 
of sensitivity, is hard (although some indication is possible), while the polar- 
ization signal is grossly affected by the overall r value. Hence, considering just 
the case of sharp reionization does not limit the significance of our conclu- 
sions. We show results obtained if the model M. has spectral index Ug — 1; 
but significant quantitative shifts occur only for true Us values unexpectedly 
different from unity. 

The SPOrt sensitivity quoted in Table2 corresponds to a pixel noise 



for an Nside = 16 pixilation; here At and Eq are the flight duration and the 
detection efficiency, respectively. Previous analysis carried out by Colombo 
& Bonometto (2003) explored the cases ap = A jiK, and up = 2 jjK. In this 
paper we report results for up — ?> as well, this value being closer to the 
SPOrt expected sensitivity. For anisotropy data, assumed to originate from 
an independent experiment, we take ctt = 2 /xK. A similar noise is typical of 
WMAP anisotropy data, once rescaled to our apertures. Cosmic and noise 
variances are taken into account by performing 1000 independent realizations 
of the model M.. 

In Figs. 19-20 we show joint confidence regions obtained in four realizations, 
chosen at random, for three different values of r in the model with cxp = 
4 fiK. as well as a realization with ap — 2 jxK. 

A plot similar to Figs. 19-20 is shown by Spergel ct al. (2003), reporting 
WMAP results. A direct comparison shows that we expect to constrain r 
better than the 1-year WMAP release, already with op = AjjK. 

Let us also stress that, on very large angular scales, information on the optical 
depth is carried almost exclusively by the polarization (and cross-correlation) 
APS. This can be seen by comparing the likelihood contours in Fig. 20, where 
the noise on T data is the same. 
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Fig. 19. Likelihood contours from the analysis of both polarization (Q and U) and 
anisotropy data, with ap = 4/liK and ctt = 2/liK, for r = 0.20 (left) and r = 0.10 
(right). We assumed a flat model with DE due to cosmological constant, J7rn = 0.3, 
h = 0.65, = 0.22, n^,. = 1. Four independent realizations of the model are shown 
(see text). More or less heavily shaded areas indicate 1 or 2-a confidence regions, 
on the Ug-T plane. 
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Fig. 20. As in Fig. 19, with r = 15 and either ap = 4;uK (left) or crp = 2 /xK (right). 
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Fig. 21. Average likelihood distribution, from the analysis of both polarization (Q 
and U) and anisotropy data, with ap = 3 jiK and ax = 2 jiK. The underlying model 
is the same as in Fig 19. 1 or 2-a confidence regions, in the n^-r plane, are shaded 
as in previous figures. 

In Fig 21 we show the likehhood contours, averaged on cosmic and noise vari- 
ances, for ap — 3 /iK. Average results are often used to describe the expected 
performance of an experiment. However, when the likelihood distribution is 
not expected to approach a Gaussian behaviour, a more detailed analysis of 
how cosmic and noise variance can affect data is in order. Accordingly, Fig. 22 
shows in which fraction of realizations various r limitations can be achieved. 
Already for ap — A //K, minimum and maximum values of r are detectable in 
90 % of realizations, at 2-a level, provided that the real r ^ 0.18. If r ~ 0.14, 
such percentage reduces to ~ 60 %. These figures should not be confused 
with those obtainable from the average likelihood. In the latter case, still for 
ap — A /iK and at the 2-a level, the minimum r value should be distinguish- 
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Fig. 22. Fractions of realizations where a lower limit on r is found, at different 
confidence levels, assuming = 2 /xK and three different values for up. 
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able from zero if the real r is just above 0.10. In this very case, there is also 
a 20 % of a-priori probability that a minimum value of r is fixed at the 3-a 
level. 

For the case ap — 3 /iK, minimum and maximum r values can be detected, in 
90% of realizations, at 2(3)-cr C.L., provided that the real r >~ 0.12(0.14). 
Such percentage reduces to ~ 60 % only if the real r is below 0.07 (0.085). 
Then, if we refer to the average likelihood, r should be lower than 0.06 not to 
be distinguished from zero. 

It is however worth pointing out that, even if r was as low as 0.03, a value 
smaller than the lower WMAP limit, we should still have a ~ 30%(20%) 
probability to distinguish it from zero at 2(3)-(T C.L. Most of these estimates 
can be deduced from Fig. 22, in which we also point out the case r = 0.17, as 
from WMAP peak probability. 

Let us remind that an improvement in r determination is certainly welcome, 
as widely discussed in Sec. 2. Large r values were the main finding of the 
first-year WMAP release and the physical imphcations of such r's concern 
the evolution of the Universe at the eve of object formation. Such results arise 
from the TE correlation at small /. Testing such correlation on data coming 
from instruments measuring T and Stokes parameters separately is, in our 
opinion, extremely important and urgent. 



5 Conclusions 

The recent results of the NASA- WMAP satellite put particular emphasis on 
the importance of CMBR-polarization investigations on large scales, able to 
either confirm or revise the estimate r = 0.17 obtained from the Ctei spectrum 
of the WMAP first-year data. Such an unexpectedly high value of r suggests 
fascinating scenarios, some of them discussed in Sect. 2.3, like a double reion- 
ization. However, due to systematics still to be fully understood, WMAP did 
not produce any polarization map yet. Thus, after their first results, there 
are more and more reasons for investigating large-scale properties of CMBR 
polarization, independently from the anisotropy. 

Previous attempts to measure the CMBR polarization had already demon- 
strated that the istantaneous sensitivity of the instruments, which essentially 
depends on the front-end noise, does not always reflect in the final sensitiv- 
ity. This was due to the presence of systematics cancelling out the advan- 
tages of long integration times. Such systematics have been widely treated for 
anisotropy experiments, but not so deeply for polarization experiments that 
led to positive detections much more recently. 
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The SPOrt experiment represents the first attempt to measure the Q and 
U Stokes parameters of the microwave sky on large angular scales as direct 
outputs of the instrument. No off-line process is needed to extract the primary 
information, e.g. the linear polarization in each pixel as defined by the 7° 
FWHM of each SPOrt antenna. The SPOrt radiometers have been designed 
to minimize every source of systematics (in polarization) following an original 
analytical approach that could lead the way to next generation instruments 
aiming at B-mode investigations of CMBR. A quite new philosophy, in fact, 
has been adopted beginning from the identification of the most critical parts 
of the instrument, passing through their requirement definition and ending up 
with custom realization of components that actually represent a new state- 
of-the-art. 

Following this idea several new devices have been reahzed, as for example the 
onboard calibrator, based on a complex of marker injectors, that allows a full 
monitoring of the instrument response. The analog correlation is performed by 
a custom passive microwave device, namely the Hybrid Phase Discriminator 
and, similarly, other critical components of the SPOrt radiometers have been 
realized reaching performances often two-three orders of magnitude better 
than equivalent off-the-shelf components. 

As far as data-reduction and analysis are concerned, a new destriping tech- 
nique has been studied to minimize the correlated noise due to residual sys- 
tematics and typically affecting the low order multipoles, where most of the 
cosmological information resides. As a consequence, the SPOrt capability to 
measure large scale CMBR polarization shall be similar to that of WMAP. 
Although it is not as high as that planned by PLANCK, SPOrt will make use 
of true polarimeters. 

In parallel to the efforts on both the analysis and the technological develop- 
ment, a relevant activity has also started for understanding and modelling the 
expected major source of polarized foreground, e.g. the Galactic synchrotron 
radiation. Existing data sets on Galactic polarized emission, unfortunately 
limited to regions around the galactic plane and frequencies below 2.7 GHz, 
have been studied in their spectral and angular properties. Since it was recog- 
nized that any simple extrapolation from few GHz to much higher frequencies 
would fail, the strategy was to build a template estimating the synchrotron 
polarized emission free from Faraday rotation. The result is a Faraday-free 
template whose computed APS can be extrapolated up to 90 GHz, at least 
for multipoles / < 20. 

In summary, SPOrt can be expected to provide a direct measurement of r, 
discriminating it from zero at the 2-a level, if the real r >~ 0.06. This as- 
sumes an excellent performance of observational apparati, leading to a pixel 
sensitivity ap — ?> fjK for a HEALPix-like pixilation with Nside — 16. Actual 
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results, of course, are subject to cosmic and noise variances. If wc consider 
that in more detail, minimum and maximum r values can be detected, in 
60 % (90 %) of realizations, at 2 a level, if the real r >- 0.075 (0.12). Finally, 
if r ~ 0.17, according to WMAP peak probabihty, not only considering the 
average hkelihood, but also in 83 % of reahzations, the detected value will be 
distinguishable from zero at the 3-a level. 

Furthermore, at the l-cx level, some separate information on the reionization 
redshift Zri and ionization rate x^, may be obtainable (see Fig. 3). In the same 
way, we can hope to put some constraints on the nature of Dark Energy, thanks 
to the effects on CMBR of different ISW patterns. The SPOrt capability to 
obtain separate information on cosmic reionization histories and DE nature, 
is however to be explored in more detail. It should never be forgotten that the 
SPOrt outputs will be completely free of any possible bias due to instrumental 
polarization-anisotropy correlation. 

Besides testing reionization and discriminating between different values of the 
Universe optical depth, which is its most appealing goal, SPOrt will provide 
almost full sky polarization maps of the Galaxy at 22 and 32 GHz. Last but 
not least, SPOrt will allow testing of new technologies presently beyond the 
state of the art, thus opening a real window upon more challenging scientific 
and civil space applications. 
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